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Abstract: It has been stated in the literature that the case of maximal mixing angle 
for v e leads to no day-night effect for solar neutrinos and an energy independent 
flux suppression of ~. While the case of maximal mixing angle and Am 2 in the 
MSW range does lead to suppression of the electron neutrinos reaching the earth 
from the sun by P$ = §, the situation is different for neutrinos that have passed 
through the earth. We make the pedagogical point that, just as with smaller mixing 
angles, the earth regenerates the \u\) state from the predominantly [1/2) state reaching 
the earth, leading to coherent interference effects. This regeneration can lead to 
a day-night effect and an energy dependence of the suppression of solar electron 
neutrinos, even for the case of maximal mixing. For large mixing angles, the energy 
dependence of the day-night asymmetry depends heavily on Am 2 . With a sufficiently 
sensitive measurement of the day-night effect, this energy dependence could be used 
to distinguish among the large mixing angle solutions of the solar neutrino problem. 
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1. Introduction 

The recent Super-Kamiokande announcement that atmospheric neutrinos are nearly 
maximally mixed has renewed much interest in the possibility that solar neutrinos 
might also be maximally mixed. In this paper we will consider only two-neutrino 
mixings, so by "maximal mixing" we are referring to the possibility that the two 
lightest mass eigenstates, |z/i) and |z/ 2 ), with eigenvalues m-i and m. 2 respectively 
(mi < m 2 ), are each equal-probability superpositions of the flavor eigenstate |z/ e ) 
(electron neutrino) and some other state \u x ), where \v x ) can be any linear combi- 
nation of \Vfj) (muon neutrino) and \v T ) (tau neutrino). Many theoretical models 
have been proposed predicting the possibility of such maximal mixing (for exam- 
ple, see [|T], 0, |3], f|, |5], |5], [JJ). In this paper we are concerned only with the MSW 
solutions to the solar neutrino problem, first proposed by Mikheyev, Smirnov, and 
Wolfenstein || |9|, |10|1 ) while the alternative possibility of nearly maximally mixed 
vacuum oscillations has been considered by other authors |TTJ. The MSW effect 



results from the neutrino interaction with matter, causing an enhancement of the 
conversion process transforming v e into v x . The MSW effect is also capable of driv- 
ing neutrinos back towards a v e state after passing through the earth. This process 
would result in a change in the v e flux between daytime and nighttime measure- 
ments, a phenomenon known as the day-night effect , or more generally zenith angle 
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dependence. Over the past decade there have been extensive studies of the day-night 
effect If 



13, 14, 15, 16, 17, 18, 19 



which have been mostly concerned with the 
small mixing angle solutions to the solar neutrino problem. Most of these studies 
have used the Mikheyev-Smirnov expression [2(J to describe the effect of the earth 
on the solar neutrinos, which we will hereafter refer to as Eq. ( |1 . 1| ) : 
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Here P$ is the probability that an electron neutrino (|^ e )) originating in the sun will 
be measured as an electron neutrino upon reaching the earth, P$e is the probability 
that an electron neutrino originating in the sun will be measured as an electron neu- 
trino after passing through the earth, P 2e is the probability that a pure \p 2 ) eigenstate 
entering the earth will be measured as an electron neutrino when it emerges, and 9y 
is the vacuum mixing angle, defined by 



l^l) = We) COsOy — \v x ) SU1 6y 

1^) = Wx) cos 9y + |z/ e ) sin6V 



(1.2a) 
1.2b) 



In the previous studies of the day-night effect several authors have claimed that 
there is no day-night effect at P5 = 1/2 [T^, [TJ|. In fact, in the formula for the day- 
night effect which is conventionally used, Eq. ( |1.1| ), the properties of the earth enter 
only through P 2e , which is explicitly multiplied by (1 — IPs)- We wish to emphasize, 
however, that the case of maximal mixing is an exception to this statement. For 
maximal mixing Eq. ( |1 . 1|) is ill-defined, since cos28y = 0, and we will show below 
that generically there is a day-night effect for this case. Nonetheless, we have no 
disagreements with either the equations or the contour plots in the aforementioned 
papers, which in fact do show non-zero day-night effects at maximal mixing. The 
purpose of this paper is to clarify the previous papers, and also to investigate more 
carefully the role of the day-night effect for maximal mixing. We will show that at 
maximal mixing P$e 7^ 1/2, implying a day-night effect and an often overlooked 
energy-dependence of the suppression of the solar neutrino flux. 

Physically, the day-night effect survives because the neutrino beam reaching the 
earth, for all MSW solutions, is predominantly |z/ 2 ). For maximal mixing this state is 
half v e and half u x , but there is a definite phase relationship, \u 2 ) = (\u e ) + \u x )) / \/2, so 
the density matrix is not proportional to the identity matrix. A coherent component 
of is regenerated as this beam traverses the earth, leading to interference with 
the incident |z/ 2 ) beam. The case is rather different from the small mixing-angle case, 
for which Eq. really does imply the absence of a day-night effect when P$ = 1/2. 
For a small mixing angle Ps equals 1/2 only when conditions in the sun drive the 
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ensemble into a density matrix proportional to the identity matrix, in which case the 
earth would have no effect. 

In the remainder of this paper we explain in more detail why maximal mixing can 
result in a day-night effect. In Section 0, we review the derivation of Eq. ( |1 . 1| ) as given 
by Mikheyev and Smirnov , and we resolve the maximal mixing ambiguity. Next 
in Section |3] we present results of numerical calculations showing the day-night effect 
at maximal mixing. Finally in the appendices, we provide greater details concerning 
the analytic and numerical calculations presented in the paper. 

2. Derivation of Equation (1.1) 

The key assumption necessary for the derivation of Eq. (|1 . 1|) is that the neutrino 
beam arriving at the earth can be treated as an incoherent mixture of the two mass 
eigenstates \i>i) and \v 2 ). That is, we assume that there is no interference between 
the v\ and v 2 components reaching the earth, or equivalently that the off-diagonal 
entries of the density matrix in the v\-v 2 basis are negligibly small. The physical 
effects which cause this incoherence are discussed in Appendix 0. In the case of 
maximal mixing, the incoherence is ensured for Am 2 > 6.5 x 1CT 9 eV 2 because of 
the energy resolution of current detectors. Other sources of incoherence include the 
separation of \vi) and \v 2 ) wave packets in transit to the earth, the averaging over 
the regions in the sun where the neutrinos were produced, and the averaging over 
the changing radius of the earth's orbit In Appendix |C| we comment on the 

regions of parameter space for which the assumption of incoherence is valid. 

Given the assumption of incoherence, we write the fractions of and \v 2 ) flux 
from the sun as ki and k2, respectively 1 . Since there is no interference, the probability 
that a solar neutrino will be measured as v e upon reaching the surface of the earth 
is given by 

Ps = h |(z/ e |z/i)| 2 + k 2 |(z/ e |z/ 2 )| 2 
= ki cos 2 Oy + k 2 sin 2 9y 

= COS 2 fly - fc 2 COs20y , (2.1) 

where we have used Eqs. ( |1.2| ) and the fact that k\ + k 2 = 1. Similarly, the probability 
that a solar neutrino will be measured as v e after passing through the earth, when 
it is no longer in an incoherent superposition of the mass eigenstates, is given by 

Pse = hP le + k 2 P 2e , (2.2) 

where P\ e (P 2e ) is the probability that a \u\) {\u 2 )) eigenstate will be measured as 
v e after traversing the earth. Finally, the unitarity of the time evolution operator 

1 For large mixing angles, sin 2 26>y > 0.5 and 5 x < Am 2 (eV) 2 < 1 x 10~ 7 , k 2 ~ 1 and 

k x ps 0. 



3 



implies that the state vectors of two neutrinos entering the earth as and |z/ 2 ) 
must remain orthonormal as they evolve through the earth and become \ux) and 
|z> 2 ), respectively. Therefore 

Ple + P 2e = \ (v e \i>l)\ 2 +\(P e \v2)\ 2 = 1 • (2.3) 

Eq. can then be obtained by using Eq. ( |2.1|) and the above equation to eliminate 
Pi e , ki, and k 2 from Eq. (|2.2|). 

From the above derivation, one can see that the singularity of Eq. (|1 . 1| ) at maxi- 
mal mixing arises when Eq. ( |2.1| ) is solved to express fo in terms of Ps- For maximal 
mixing Ps = 1/2 for any value of fo, so fo cannot be expressed in terms of P$. The 
ambiguity disappears, however, if one leaves fo in the answer, so Eq. (|2.3| ) can be 
used to rewrite Eq. (|2.2|) as 

P SE = \ + 2 (fa - i) (P 2t - i) . (2.4) 

Thus, Pse = 1/2 only if k 2 = 1/2 or P 2e = 1/2. For the MSW solutions at maximal 
mixing one has fowl, and there is no reason to expect P 2e = 1/2. Generically 
Pse 7^ 1/2 for the case of maximal mixing. 

3. The day-night effect at maximal mixing 

Using the evolution equations derived in Appendix [A] and the procedures described 
in Appendix [B|, we have calculated a variety of properties concerning the day-night 
effect for maximal mixing angle. The calculation parameters are chosen for those of 
the Super-Kamiokande detector. 

Figure [I] shows the evolution of P{v e — > u e ), the probability that a solar neu- 
trino will be measured the beam of neutrinos traverses a path through the 
center of the earth. Notice that after traversing the earth the ensemble of neutri- 
nos is no longer in a steady state, but instead P{y e — > v e ) continues to oscillate 
in the vacuum. From the perspective of the mass eigenstates, the neutrinos under 
consideration arrive at the earth roughly in a \v 2 ) state. Upon reaching the earth, 
the step-function-like changes in the electron density profile (see Fig. [|) cause non- 
adiabatic evolution, regenerating the \v\) state and leading to interference effects. 
In the regions of parameter space where the day-night effect is maximal because the 
oscillation length of these interference terms coincides with the length of the slabs of 
near constant density composing the earth, the resulting buildup of u e flux has been 
called oscillation length resonance |22| , |23| , |24| . 

In Fig. |2| we present a contour plot calculated from the density matrix that 
exemplifies the non-zero nature of the day-night effect at maximal mixing. On the 
left is a three-dimensional surface where the height of the surface is the day-night 
asymmetry. Notice that the exposed edge is calculated at maximal mixing and is 
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Evolution of Maximally Mixed Solar Neutrino Ensemble Crossing the Earth 
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Figure 1: The evolution of P Ve -> Ve as the ensemble of neutrinos propagates across the 
center of the earth. The neutrinos enter the earth as an incoherent mixture of the energy 
eigenstates v\ and V2 which is almost completely v^. This plot shown is for Am 2 = 
1.3 x 10~ 5 eV 2 and a neutrino energy E = 6.5 MeV. 




sin 2 (20 v ) 




sin 2 (20 v ) 

Figure 2: The day-night asymmetry (A^-n = (N — D)/(N + D)) as a function of mixing 
parameters calculated using the density matrix. On the left is a three dimensional surface 
where the height of the surface is the day-night asymmetry. Notice that the exposed edge 
is calculated at maximal mixing and is clearly non-zero. On the right is a contour plot 
showing the lines of constant day-night asymmetry. 



clearly non-zero. On the right is a contour plot showing the lines of constant day- 
night asymmetry, a plot which is identical to those produced in other references. 

We now explain how the inclusion of the day-night effect at maximal mixing 
resolves a certain confusion that has arisen in the past because of its neglect. As a 
result of the non-zero day-night effect, there exists an energy dependence at maximal 



5 



Energy Dependent Flux Supression at SuperKamiokande 
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Figure 3: The predicted flux suppression as a function of energy. Notice that the predicted 
overall flux suppression is not 1/2, due to day-night effects, even though the mixing angle 
is maximal. The plot is for Am 2 = 1.0 x 1CP 5 eV 2 which is near the border of the region 
excluded by the small day-night effect (Ad- n ) measured at Super-Kamiokande. 



mixing, as can be seen in Fig. ^. If one assumes that the flux suppression at maximal 
mixing has no energy dependence, as was done in Ref. ||25| , then there is an apparent 
discrepancy between two sections of Ref. . Sec. IV-D excludes the possibility of 
energy-independent oscillation into active (as opposed to sterile) neutrinos at the 
99.8% confidence level, while Fig. 2 shows some regions of the maximal-mixing- 
angle parameter space not excluded at the 99% confidence level. Ref. [^5] has tried 
to resolve this discrepancy without including the day-night effect, concluding that 
maximal mixing is excluded at the 99.6% confidence level. The actual resolution to 
this apparent discrepancy is that Fig. 2 of Ref. |26| includes the energy dependence 
induced by the day-night effect at maximal mixing, while Sec. IV-D discusses the 
case of energy-independent flux suppression and does not apply to maximal mixing. 
The correct conclusion is that of Fig. 2, which shows that maximal mixing is not 
excluded at the 99% confidence level. 

Whether or not the day-night effect is included, maximal mixing is not a very 
good fit to the experimental data from the three neutrino experiments (chlorine, 
gallium, and water) J26j. However, maximal mixing does fit well if the chlorine data 
is excluded on the suspicion of some systematic error [^Tfl. Ref. || has argued that if 
the 8 B flux is about 17% lower than the standard solar model (BP98) |28fl , then a bi- 
maximal mixing scenario becomes a tenable solution to the solar neutrino problem. 
The MSW mechanism described here is applicable for Am 2 > 6.5 x 10~ 9 eV 2 . In the 
bi-maximal mixing scenario that we consider the upper bound on Am 2 is set by the 



CHOOZ data constraining Am 2 < 0.9 x lfr 3 eV^ p 
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Figure 4: The day-night asymmetry (Ad- n = (N — D)/(N + D)) as a function of recoil 
electron energy at Super-Kamiokande. Both plots are at maximal mixing angle, with Am 2 
at the upper and lower borders of the region disfavored by the smallness of the day-night 
effect observed at Super-Kamiokande. The rising line is for Am 2 = 2 x 1CP 5 eV 2 , and the 
descending line is for Am 2 = 3 x 10 -7 eV 2 . 



When detailed studies of the day-night effect are completed, the energy (and 
zenith angle) dependence will be valuable additional information. To the best of 
our knowledge, the Super-Kamiokande Collaboration has not published their day- 
night asymmetry as a function of recoil electron energy. Past studies of the day- 
night effect have noted the energy dependence of the day-night asymmetry [[H], |T8[ |. 
While for small mixing angles |A^_ n | < 0.02 without a clear energy dependence 
[18] , for large mixing angles the Ad_ n energy dependence can be significant and 
informative. Fig. ^ shows the theoretical predictions of the day-night asymmetry in 
the electron recoil spectrum at Super-Kamiokande for two cases of maximal mixing: 
Am 2 = 2 x 10™ 5 eV 2 and Am 2 = 3 x 10~ 7 eV 2 . Note that the two curves have 
opposite slopes. 

The approximate shape of the graph of Ad- n vs. recoil electron energy can be 
understood from Fig. |2], using the fact that Fig. |2] is dominated by the peak of the 8 B 
neutrino spectrum at about 6.5 MeV. It is shown in Appendix [A] that the neutrino 
evolution equations (Eqs. ( |A.5|) -( [A!8|) ) depend on Am 2 and the neutrino energy (or 
momentum) E only through the combination Am 2 / E. Thus, Fig. |2| shows that for 
any value of sin 2 26>y, A d _ n has a maximum at Am 2 / E « 2.5 x 10~ 6 eV 2 /(6.5 MeV). 
When E is varied at fixed Am 2 , Ad~ n will have a peak at 



E 



Am 2 



2.5 x 10" 6 eV z 



x 6.5 MeV 



(3.1) 



So for Am 2 = 2 x 10 5 eV 2 the peak lies far to the right of the scale in Fig. f|, so the 
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curve slopes upward. For Am 2 = 3 x 10~ 7 eV 2 the peak lies far to the left, and the 
curve slopes downward. 

Fig. ^ shows that the peak in the graph of A^-n vs. Am 2 is higher at large 
mixing angles (sin 2 2#y « 0.7) than it is at maximal mixing, so the same will be 
true for the energy dependence of the day-night effect. For sin 2 2#y = 0.63 and 
Am 2 = 1.3 x 10~ 5 eV 2 , for example, the slope of the graph of Ad- n vs. recoil electron 
energy is about twice the magnitude of the slopes shown in Fig. £|. Thus, the day- 
night asymmetry as a function of recoil electron energy could be a strong indicator of 
Am 2 if the solar neutrinos have a large or maximal mixing angle in the MSW range 
of parameters. 

4. Conclusions 

We have disproved the assumption that Pg = 1/2 always implies P$e — 1/2. We have 
also shown that neutrinos with a maximal mixing angle can have a day-night effect 
and that they do not always result in a uniform energy-independent flux suppression 
of 1/2. Because the issues that we have attempted to clarify concern mainly the 
words that have been used to describe correct equations (which were generally used 
numerically), there are no changes to most constraints presented in other references. 
The only corrections apply to fits of energy-independent suppressions; that is the 
fits no longer apply to the exclusion of some regions of maximally mixed neutrinos. 
Finally, we have noted that the energy dependence of the day-night effect can be a 
strong discriminator between various solutions of the solar neutrino problem. 
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A. Derivation of MSW equations 



In this appendix we present a derivation of the MSW effect. First we derive the 
MSW equations of motion for an individual neutrino. We then find the energy 
eigenstates of the system and use them to find the wave function amplitudes for 
electron neutrinos produced in the sun and evolved into the vacuum. To describe 
the ensemble of neutrinos we introduce the density matrix. After averaging out the 
rapid oscillations we find a steady state solution to the density matrix equations of 
motion. We average this solution over the regions of neutrino production. 

We begin by finding the MSW equations of motion for an individual neutrino. 
The coupling describing the interaction between electron neutrinos and electrons is 

H mt = V2G F N e , (A.l) 

where N e is the number density of electrons. This contribution to the interaction 
Hamiltonian is added to the Schrodinger equation written in the flavor basis. We 
assume that \u e ) can be written as a superposition of only two mass eigenstates, \ui) 
and \u 2 ). We let \v x ) denote the orthogonal linear combination of and \u 2 ), which 
might be any superposition of |z/^) (muon neutrino) and \u T ) (tau neutrino). The 
transformation between the v\-v 2 and u e -u x bases is then given by 

(C V1 \ = ( cos6 v - sin fly \ / C v \ 

\C U J \sm6 v cos9 v J \C U J ' 1 ' ' 

where the variable 9y is the vacuum mixing angle, and C„ = for v = vi, u 2 , u e 

or u x . This equation can be written compactly by introducing the index notation 

C Vi = Uj f C Uf , (A.3) 

where the repeated index / is summed over v e and v x , and % is summed over the 
mass eigenstates. The Schrodinger equation for this system is: 





where we have expanded the energy in the ultra-relativistic limit so that E — p+ 
We now substitute U into the Schrodinger equation, obtaining 





where 



(A.5) 



A = ^-(m 2 2 - ml) (A.6) 
A=^-sm26 v (A.7) 
B = ^G F N e - ^cos2fl y , (A.8) 
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and where we have dropped the term p+ {jn " L ~^™ 2) + ^GpN e which is proportional to 
the identity, because terms proportional to the identity cannot contribute to mixing. 
The eigenvalues are ±X(N e ), where X(N e ) = \J A 2 + B 2 , and the eigenvectors 

W) ' and " + = (v#J' (A ' 9) 

Since these eigenvectors form the matrix that will diagonalize the interaction matrix 
in the presence of matter, it is useful to parameterize them by a matter mixing angle 
M (N e ): 

cos 9 M = J , and sm6 M = J ^ , (A.IO) 



or equivalent ly 



or 



Defining the matrix 



Xcos26 M = -B, (A. 11) 

A sin 29m = A. (A. 12) 



^«)=( C °^' (A.13) 



sin^M cosQm 
the Hamiltonian can be diagonalized as 




U\e M ) = ( J . (A.14) 



We maintain the notation introduced in Eq. ( |A.3| ) so that C Vi (9m) = UL(0M)C Vf in 



or out of matter, where C^Om) = (^1^) denotes the amplitude for the overlap of 
the neutrino state with instantaneous mass eigenstates 

To describe the evolution of the neutrinos as they travel to the earth from their 
creation point in the sun, it is useful to develop the adiabatic approximation, in which 
one assumes that the density changes imperceptibly within an oscillation length. 
Remembering that U, 6m, and A are all functions of the local electron density A^ e , 
and hence functions of time, we write the Schrodinger equation in the basis v^Om) 
of the instantaneous mass eigenstates: 

Mc:)-(~o °)(ci;) + («(a;) (A15) 



-A 









I) 


r 


-A 




IC 





I) 


c, 



V 2 





(A.16) 



The adiabatic approximation is the assumption that the off-diagonal terms 8 x 6m can 
be neglected, in which case the equation is easily integrated: 



C^tf)\^(e + ^ \(C Ul (to) 



(A.17) 
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where 



to 



\(t)dt . 



(A.18) 



Because the adiabatic states form a complete basis, we can always write the exact 
solution as a superposition of the two adiabatic states. This final superposition is 
expressed by two unknown variables, a x and a 2 where |ai| 2 + \a 2 \ 2 = 1. The | a-2 1 2 
parameter represents the probability of a non-adiabatic transition, which is most 
likely to happen when the neutrinos cross resonance, the density at which B = 0, 
when the two eigenvalues become nearly equal. Likewise | «i | 2 = 1 would represent 
adiabatic evolution. Given any initial state Vf(to) in the flavor basis, the final state 
can be written in the general form: 



c V2 (t f ) 











m )tf(e M (t )) v f (to) . (A.19) 



For an electron neutrino originating in a medium of mixing angle 9m, the above 
equation implies that the final state in the vacuum is given by 



'c n (t f y 

C V2 (t f ) 




ai cos 9 M e +l ^ + a 2 sin 9 M e ^ 
-a 2 cos 9 M e +i ^ + a\ sin 9 M e~^ 



(A.20) 



We now go on to talk about the ensemble of neutrinos reaching the earth. To 
describe a quantum mechanical ensemble of neutrinos, it is useful to introduce the 
density matrix 

p = J2fiWi)(ui\ , (A.21) 

% 

where fi denotes the probability that the particle is in the quantum state |z/j). The 
density matrix corresponding to a single neutrino as described by Eq. (|A.20| ) is 
therefore given by 

'\A X \ 2 AiA* 2 \ 



P 



A\A 2 \A 2 \ 2 



where 



A X A\ 



\A,? 



1 + cos26» A/ (l - 2\a 2 \ 



+ 



sin 29 



2 -2i<t>{t f ) 



a e 



cos 29 



M 



2\a 2 \ 



a x a\ sin 29 M e 2i<i,{t f ) + c.c 

— aia 2 cos 29m 

a x a 2 sin 29 M e 2i ^ ) + c.c 



(A.22) 

(A.23) 
(A.24) 
(A.25) 



In Appendix we explain why this process allows us to eliminate the terms that have 
rapidly oscillating phases. In particular, the phase angle </>(£/) and the phases of the 
complex numbers a\ and a 2 are all rapidly varying functions of the neutrino energy, 
the location in the sun where the neutrino is produced, and the precise time of day 
and year at which the neutrino is observed. The density matrix which describes the 
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ensemble of observed neutrinos is constructed by averaging over these quantities, so 
any quantity with a rapidly oscillating phase will average to zero. This is equivalent 
to the statement that the v\ and components arriving at the earth are incoherent, 
so we average over their phases. The matrix elements of the phase-averaged density 
matrix are given by 



(|Ax| 2 ) = - [l + cos2fl A/ (l - 2|a 2 | 2 )] (A.26) 
(AxAfi = (A.27) 
(\M 2 ) = \ [l-cos2^ / (l-2|a 2 | 2 )] . (A.28) 



The term \ai\ = -fjump is the probability of crossing from one adiabatic state to 
the other during the time evolution of these operators. An approximate expression 
for Pj Ump can be found by using a linear approximation for the density profile at 
p0| , yielding 



resonance 



\ Apcos(20 v )N'(x rca ) J 

Here N(x res ) is the density at the point where the neutrino crosses resonance, and 
N'(x Tes ) is the first derivative of the density at resonance. More accurate approxima- 



tions to Pjump and the details of their derivation can be found in Refs. [[Jl], [32| and 
the references therein. 

The density matrix corresponding the ensemble of observed neutrinos must be 
obtained by averaging over the production sites in the sun. While we have already 
made use of this fact in dropping all terms with rapidly oscillating phases, we must 
still average the slowly varying terms which remain. Letting 8 B(r) denote the nor- 
malized probability distribution for production at a distance r from the center of the 

(A.30) 



sun, one finds 




where 

Co = lf dr s B(r) cos(29 M (r)) (1 - 2P jump ) . (A.31) 

Note that the diagonal entries of p are just the fractions k\ and /c 2 of v\ and v-i flux 
from the sun, respectively, as defined in Sec. 0. Therefore 

h = ~ + C , k 2 = ~ - C . (A.32) 
Finally, we transform to the v e -v x basis, so 

= u(e v ) P u\e v ) = ( pee pex ) . (A.33) 
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One then finds that the probability of observing a neutrino reaching the surface of 
the earth as an electron neutrino is given by 

Ps = Pee = i + Co COS 20y . (A.34) 

The off-diagonal matrix element is given by 

p xe = — C sin 26V • (A. 35) 

Our numerical simulations were all performed by integrating Eq. ( |A.5|) to solve 

for P 2 ei and also by integrating the density matrix equations of motion. The evolution 
of the density matrix is given by 

ihdtp = -[p, H}. (A.36) 



Using Eq. ( |A.36| ) with the Hamiltonian in the flavor basis, we find that our new 



equations of motion are 

idtpee = A(p xe - p* xe ) (A.37) 

id tPxe = 2{A Pee - B Pxe ) - A , (A.38) 

where A and B are defined in Eqs. (|A.7|) and (|A.8|) . This allows us to perform 
calculations using the complete mixed ensemble. The expressions given in Eqs. ( |A.34j ) 



and (|A.35|) form a steady state solution of the density matrix equations of motion 



in the vacuum. For typical (Am 2 m 1 x 1CT 5 eV 2 and E = 8 MeV) maximally 
mixed Boron-8 ( 8 -B) neutrinos, P 2e oscillates and can take any value between and 
1. Co ~ — \, which means that k 2 ~ 1, and Eq. (|2.4| ) reduces to Pse = Pie- Thus 
Pse exhibits oscillatory behavior, and is in no way constrained to be 1/2 at maximal 
mixing. 

B. Calculation methodology for the day-night effect 



First we calculated Pg using Eq. ( |A.34| ) for the spectrum of Am 2 /p at various mixing 
angles. For a given Am 2 /p we averaged P$ over the regions of S B neutrino production 



in the sun, provided by Ref. [Q. Using Pg to describe the neutrinos that arrive at 
the earth, we then performed the evolution through the earth with the density matrix 
equations of motion. The initial conditions for the density matrix are given by 

Pee = P S , (B.l) 

and 

p xe = --(2P s -l)taa20 v . (B.2) 

At maximal mixing we assume that p xe = — | cos(^26 l j\/(to)) sin 26V ~ § which is the 
adiabatic result. This assumption is justified because in the regions of parameter 
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Figure 5: The Preliminary Reference Earth Model (PREM) electron density (N e ) profile 
of the earth. N e is shown in units of Avogadro's number of electrons per cm 3 . 



space under consideration near maximal mixing, Pj um p ~ 0. It follows that in these 
same regions of parameter space the evolution remains adiabatic in the limit where 
Qy = 7r/4. We use the earth density profile given in the Preliminary Reference 
Earth Model (PREM) []33| (see Fig|5|). To convert from the mass density to electron 
number density we use the charge to nucleon ratio Z/A = 0.497 for the mantle 
and Z/A = 0.467 for the core. The numerical calculations were performed using 
a fourth order Runge-Kutta integration programmed in C++. We propagated the 
neutrinos through the earth for 90 zenith angles, a, evenly spaced between 90 and 
180 degrees. We calculate the anticipated electron flux as a function of zenith angle 
and energy, denoted Pse{&, E u ). The calculation parameters are chosen for those of 
the Super-Kamiokande detector. The normalized 8 B neutrino spectrum, <&(E U ), and 



solar electron densities, N e , are also obtained from data-files provided by Ref. [33 



Effective neutrino cross sections are available which take into account the electron 
recoil cross section with radiative corrections, the energy resolution, and the trigger 
efficiency |S3[. We used these more accurate cross sections for the overall day- 
night effect plotted in Fig. ||. Because these effective cross sections already include 
the integration over detected electron recoil energy, to calculate the recoil electron 
spectrum we used the differential neutrino-electron scattering cross sections given 



in Ref. [36]. Using these data files and numerical results the cross section for the 
scattering of solar neutrinos of energy E v with electrons to produce a recoil electron 
of energy T' at the zenith angle a is given by 

d -^(T\E^a) = P SE (a,E^(T',E u ) + [1 - P SE (a,E v )} ^(T',E U ). (B.3) 



Since muon and tau neutrinos have the same neutral current interactions, we can use 
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the Up cross section for the v x . The analysis of the recoil electron spectra is explained 
in Refs. f35| . The actual flux at recoil energy, T, is 

g(a,T) = j o dE v *{E V ) jf rfT' i?(T,T') ^(1*^0) (B.4) 
where the energy resolution of the detector is incorporated through 



The energy resolution, A^/, around the true electron energy T' for Super-Kamiokande 
is given by 



At' = (1.6 MeV) y/T'/(10 MeV). (B.6) 

To calculate the average day-night effect over one year, we weight the flux by the 
zenith angle exposure function Y(a) explained in Appendix [B]. The daytime mea- 
sured flux at a given measured electron recoil energy, T, is given by 



/■90 

D{T)= da g(a,T)Y(a), (B.7) 
Jo 

and for nighttime is 

r 180 

N(T)= da g(a,T)Y(a). (B.8) 

J90 

The day-night asymmetry as a function of recoil electron energy plotted in Fig. |] is 
given by 

A m _ N(T) - D(T) 

A d - n {T) - N(T) + D{T y (B.9) 
The final day-night asymmetry plotted in Fig. |2| is given by 

ir Me vdT (N{T)-D(T)) 

A d _ n = y— — f (B.10) 

Is°°u e vdT(N(T) + D(T)) 



where 5 MeV is the minimum energy detected at Super-Kamiokande. Verification of 
the accuracy of the computer code has been accomplished with the help of |37| , and 
by comparing our simulations to plots and data available in the literature. 



C. Validity of the steady state approximation 

Most of the work in the past decade on the MSW effect has assumed that the ensemble 
of neutrinos reach the earth in a steady state solution of the density matrix (i.e., in 
an incoherent mixture of the mass eigenstates v\ and u 2 ). There are several reasons 
that the neutrinos reach the earth in a steady state: (a) The separation of the 1^) 
and ^2) wave-packets while propagating from the sun to the earth exceeds the size of 
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the individual wave packets, eliminating the interference effects, (b) The eccentricity 
of the earth's orbit results in a daily change of the earth-sun radius larger than the 
vacuum oscillation length of the neutrinos, (c) The neutrinos are produced in a 
region much larger than their local oscillation length, (d) The energy resolution of 
the current detectors coupled with the earth-sun radius perform an average. We now 
proceed to map out parameter space justifying where the steady state approximation 
is valid. First we consider the separation of the two eigenstates during transit to the 
earth. This results in system that is an incoherent superposition of and l^). 
The width of the wave-packets, a x , is given by Ref. [08| : 



a x « 0.9 x 1(T 7 cm. (C.l) 

This results in a coherence length given by: 

r IE 2 

L C oh = 2V2a x - — -. (C.2) 
Am A 

We lose coherence between the mass eigenstates if L co h < 1 AU = 1.5 x 10 13 cm. 
If we require that the incoherence condition apply up to 14 MeV to include all 8 B 
neutrinos, we find that for all of sin 2 26V where Am 2 > 6.63 x 1CT 6 eV 2 the wave- 
packets have separated upon reaching earth. This corresponds to the region labeled 
(a) in Fig. |]. Because there is a continuous beam of neutrinos arriving from the sun, 
we can ignore the fact that the lighter mass eigenstate arrives first, and simply drop 
terms that rapidly oscillate due to the lack of interference between the two states. 

In the previous case the interference effects vanish because of a loss of coherence 
between the mass eigenstates for a neutrino produced at a specific place and time. In 
the remaining topics the interference effects vanish due to averaging over the ensemble 
of neutrinos which reach the detector. Regardless of the source of the incoherence, 
any two systems with the same density matrix behave identically in single counting 



experiments [39]. 



Next, we analyze the effect of the eccentricity of the earth's orbit . We are 
interested in day-night effects; therefore, if the earth-sun radius changes more than an 
oscillation length during one day, this will result in washing out any phase dependence 
in the results measured over a period of one year. Between perihelion and aphelion 
the earth-sun radius changes 2e(l AU) = 5.1 x 10 11 cm, where e = 0.017 is the earth 
eccentricity. The earth-sun radius changes by this quantity once each 180 days giving 
an average daily change in radius of 2.83 x 10 9 cm. This ensures our incoherent phase 
for Am 2 > 1.2 x 10" 6 eV 2 . This region is denoted by everything above the line 
marked (b) in Fig. |6|. 

Third, we study the impact of where the neutrinos were produced. If the neutrino 
region of production is greater than the local oscillation length of the neutrinos, then 
neutrinos of all possible phases exist in the ensemble. For a continuous beam of 
neutrinos, this also results in dropping the rapidly oscillating terms. The condition 
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is satisfied for the entire parameter space under consideration 0.001 < sin 2 2#y < 1 

"I "1 2 9 2 

and 1 x 10 eV < Am < 1 x 10 eV . However, one must be careful in 
making this statement. Although the region of production may be greater than the 
neutrino oscillation length in the sun, the neutrinos could undergo a non-adiabatic 
transition bringing a specific phase into dominance. This is the case for vacuum 
oscillations (Am 2 m 4 x 10~ 10 eV 2 ). The 8 B neutrinos are produced mostly at 
Rsb = 0.046 R sun = 3.2 x 10 9 cm. The vacuum oscillation length is on the order of 
1 AU. However the oscillations length near the solar core where these neutrinos are 
produced is about 1.8 x 10 7 cm <C Rs B . Although the neutrinos are produced in a 
region larger than their oscillation length, they acquire roughly the same phase in 
the process of leaving the sun. This occurs because the density change upon leaving 
the sun occurs more rapidly than the oscillation length of the neutrinos, violating 
the condition of adiabaticity. To express this quantitatively we estimate that if 
-Pjump < 0.1 for 14 MeV neutrinos that the initial randomly distributed oscillation 
phases at the time of production will persist as the neutrinos leave the sun and enter 
the vacuum. This leads to a steady state solution applicable in the parameter space 
above the diagonal line labeled (c) shown in Fig. || 

Last, we study the impact of the energy resolution on our ability to discriminate 
phases. Assuming perfect coherence between the two mass eigenstates the phase 
upon reaching the earth is given by 

*=^£M. ( C3, 
Apnc 

Our uncertainty in energy impacts our uncertainty in phase through error propaga- 
tion: 

Am 2 (lAin 

5P = AfKc 5P - (C ' 4) 



dp 



If our uncertainty in our phase is greater than 2n we are again justified in treating 
our ensemble as a steady state. Using conservative figures for energy (p = 14 MeV), 
and the energy resolution (Sp ~ 1 MeV) p5 |, we find that for Am 2 > 6.5 x 10~ 9 eV 2 



we are justified in the steady state approximation. This corresponds to parameter 
space above the line labeled (d) in Fig. ||. Recently Ref. |2lJ] also reached the same 
conclusions outlined in this appendix. 



D. The zenith distribution function 

The zenith angle distribution function gives the fraction of the time that the sun is 
at a given zenith angle. The function is calculated by numerically simulating the 
orbit of the earth around the sun. We begin by writing the vector towards the zenith 
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Regions Satisfying Steady State Density Matrix 



1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 




a - 




— h. : 








d J 







io- 3 
J io~ 4 



K > 

5 CD 



0.001 



0.01 0.1 

sin 2 (20 y ) 



IO" 6 

J io- 7 
io- 8 

IO" 9 
10 -io 



Figure 6: The regions satisfying the conditions for steady state density matrix. Above 
the line (a) is in steady state because of wave packet separation. Above the line (b) can be 
treated as steady state because of the eccentricity of the earth's orbit. Above the diagonal 
line (c) is in steady state because the region producing the neutrinos is much larger than an 
oscillation length, and this phase averaging survives until the neutrinos reach the vacuum. 
Above line (d) is in steady steady state because of the energy resolution of our detectors. 
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Figure 7: The zenith distribution function at Super-Kamiokande. 



of the detector in coordinates for which the earth's orbit lies in the x-y plane: 



/l 





cos 5 — sin 5 
sin 5 cos 5 



'sin(90° - L) cos. 
sin(90° - L) sin< 
cos(90° - L) 



(D.l) 



where the north latitude is given by L, <fr gives the time of day in radians, and 
5 = 23.439° is the earth's declination HO. Because we are averaging over a one year 
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time period we can arbitrarily choose the initial time of year, and the initial time of 
day. The vector pointing from the earth towards the sun is 




(D.2) 



where D is the day of the year in radians. From here we can find the local zenith angle 
from the dot product f s -f z = cos a. To numerically calculate the zenith function 
distribution we divided a into 360 bins evenly spaced between and tc. Now we run 
< D < 2n and < <p < 2n over 1000 steps in D and 1000 steps in </> and count how 
much relative time a spends over each bin. We generate the zenith angle distribution 
function for Super-Kamiokande which sits in Gifu Prefecture, Japan at 36.43° north 
latitude [|19j. This produces the undistorted zenith function distribution seen in 
Fig. [7[ One can also obtain this function as a data file from which includes small 
corrections for the eccentricity of the earth's orbit and the wobble of the earth's 
declination. To maximize accuracy we performed our calculations using this data 
file. 
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